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We investigate the spectral flow of the integrable non- 
hermitian Heisenberg spin chain under boundary conditions 
with complex twist angle. It is shown that the period of 
the spectral flow is 4n up to a certain critical imaginary 
twist, beyond which the period jumps successively to higher 
values. We argue that this phenomenon caused by non- 
hermitian properties of the system is closely related to the 
metal-insulator transition caused by non-hermitian hoppings 
for the one-dimensional insulator. 
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I. INTRODUCTION 

Boundary conditions have been playing an important 
role in statistical mechanics. Imagine for instance a 
one-dimensional (ID) many particle system with twisted 
boundary conditions where a particle gets additional 
phase factor e J * when it moves once along the system 
with a ring geometry and return to the previous posi- 
tion. For integrable systems, this general phase factor 
is simply incorporated into the Bethe ansatz equations, 
and their characteristic properties have been intensively 
studied exactly 0-0]. Physically, the twisted boundary 
condition is equivalent to introducing a magnetic flux 
threaded in the ring system ^ . In proportion as the flux 
increases, the spectrum of the system flows in a complex 
way, experiencing many level crossings and level repul- 
sions. We can get detailed information for the physical 
properties from the spectral flow thus obtained. For ex- 
ample, a local behavior of the spectral flow near $ ~ 
directly gives the transport coefficient such as the Drudc 
weight, which is a sensitive probe of the metal-insulator 
transitio n ]u J7[ . From a more global structure of the spec- 
tral flow |8|H10|, especially paying attention to its period 
for the ground state, one can also determine whether 
or not a bound state form s |TT|, r^]. Berry's phase was 
calculated in this process [Op], which may carry the 
information for statistics of elementary excitations. It 
was indeed shown that the period of the spectral 



flow is directly related to Haldane's fractional exclusion 
statistics [H. Also, the twisted boundary condition has 



attracted other current interests p9|j20| ]. 

Normal twisted boundary conditions can be further 
generalized to those with complex twist angles, $ — > 
$ + ||]. They make the Hamiltonian non-hermitian, 
which is the most remarkable difference from normal 
twisted boundary conditions. The spectral flow is now 
defined on the complex plane and shows a rich structure. 
An interesting point to be emphasized is that this prob- 
lem has the close relationship to recent topics studied 
actively in the context of non-hermitian systems. In par- 
ticular, it was claimed by Hatano and Nelson ^ij that 
introducing non-hermitian hopping terms into the Hamil- 
tonian brings about the localization-delocalization tran- 
sition even for ID random systems. The key point is that 
eigenvalues become complex if we increase the imaginary 
twist angle, which can provide a criterion whether the 
delocalization happens. Many interesting problems on 
the non-hermitian Hamiltonian systems still remain to 
be explored. In order to study such problems, integrable 
non-hermitian models are quite useful. 

In this paper, we investigate a ID non-hermitian sys- 
tem with twisted boundary conditions, and discuss how 
the spectral flow changes its character if the imaginary 
twist is introduced. As a typical example, we examine 
a version of the asymmetric Heisenberg spin chain and 
the related vertex models studied by p2|-|30|]. It was ar- 
gued that asymmetry causes interesting phenomena. For 
example, by introducing non-hermitian XY components, 
the massive XXZ spin chain with Ising-type anisotropy is 
driven to the massless phase, which is regarded as a kind 
of the metal-insulator transition (2^] . The model we will 
study in this paper is the isotropic Heisenberg spin chain, 
which is massless even for zero imaginary twist. However, 
if we consider a finite-size system, there exist many small 
gaps due to finite-size effects. It is these gaps that deter- 
mine the period of the ground state in the global spectral 
flow. We show that introducing the non-hermitian imag- 
inary twist may collapse these gaps, changing the period 
to a larger one. It is argued that this phenomenon has 
essentially the same origin as for the bulk metal-insulator 
transition in ID non-hermitian many-body systems: In 
both cases, the transition is triggered by complex eigen- 
values due to non-hermitian hopping terms. 

This paper is organized as follows. In the next section, 
we introduce the model and the Bethe ansatz equations, 
and briefly summarize the characteristics in the spectral 
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flow of the ordinary symmetric (hermitian) model. In sec- 
tion III, we discuss the spectral flow of the asymmetric 
model with complex twist angles, by numerically evalu- 
ating the Bethe-ansatz spectrum for small systems. We 
show that there is indeed the critical value of complex 
twist, at which the period of the spectral flow suddenly 
changes. In IV, we estimate this critical value for the 
transition for large but finite-size systems. For this pur- 
pose, we propose an approximate treatment to obtain the 
finite-size corrections rather accurately. Final section V 
is devoted to summary and discussions. 



II. MODEL AND BASIC PROPERTIES OF THE 
SPECTRAL FLOW 

In this section, we introduce the model and summarize 
some basic properties of the spectral flow for the ground 
state of the hermitian (symmetric) model. 



A. Model and Bethe ansatz equations 

The model we consider is the well-known isotropic 
Heisenberg spin chain 



but with the twisted boundary condition, 



°j+N — e J jj °j+N — °j ■ 



(2.1) 



(2.2) 



The spectral flow has been discussed so far for real 8. 
What we are now interested in is how the ground state 
behaves for the complex twist angle, 



e = $ + 



(2.3) 



which makes the Hamiltonian (2.1) non-hermitian. In 
the case of purely imaginary twist angle, the model is 
the same as the asymmetric Heisenberg chain but with 
the replacement $ — * Nfy. The asymmetric XXZ Heisen- 
berg chain has been studied so far in the massive regime, 
while our isotropic model is in the massless regime for 
which the non-hermitian effects disappear in the ther- 
modynamic limit. We shall in turn treat the model for 
a finite-size system, and discuss how the non-hermitian 
property affects its spectral flow by examining the X/N 
corrections. 

We begin by summarizing the symmetry properties of 
the model. Define the eigenvalue problem as 



H(Q)\^(Q)) R = E R {G)\i/>(e)) R 

L (V(e)|#(e) = L {i>(e)\E L (e) 



(2.4) 



Because of the symmetries H*(G) = H(-0*), W(0) = 

H(Q*), and UH{@)W = H(-Q) with U = e~™^ s s 
we have 

m-e*)) R = \r(e)) R , E R (-e*) = E* R (e), 

\m*))R = U\r(e)) R , E R (Q*) = E* R (G), 
MQ))l = U\1>{-Q')) r , E L (e)=E R (Q). (2.5) 

As shown by Yung and Batchelor ||, the model is 
solved by the Bethe ansatz method, 

M 

Npfa) + - A fc) = 27rI J + e ' j = 1, ■ • • , M, 

(2.6) 

where the momentum p{X) and the two-body phase shift 
9(X) as functions of the rapidity A are given by 



P(A) 
0(A) =*ln 



i In 



X-i/2 
~ \ + i/2j ' 



A + i 



(2.7) 



Here the branch cut lies on (— ioo, —i/2) and (i/2, ioo) for 
p(X) and on (— ioo, — i) and (i,ioo) for 6(X). The total 
energy of the system is given by the summation of single 
particle energies, 



M 



4Af + 1* 



(2.8) 



In the rest of the paper, we concentrate on the spectral 
flow of the ground state, so that we set M — N/2 with 
N = even, and 



N 1 

T ~ 2 



N 1 

T ~ 2 



(2.9) 



B. Spectral flow for the case $ = 

Before we proceed to the non-hermitian model, let us 
briefly summarize the spectral flow of the model with an 
ordinary hermitian condition = 0. This case was in- 
vestigated in detail by P,p|,p|, p^|Jl^ 1 , who claimed that the 
ground state becomes the first excited state at $ = 2tt 
and returns to the ground state again at $ = 4-7T, namely, 
the period of the ground state is Att as a function of $. 
This feature can be observed indeed in Fig. [l], which is 
obtained by the exact diagonalization for the S z = sec- 
tor of the TV = 6 system. The period 4ir is caused by 
the small gaps made by level repulsions due to irrelevant 
backward interactions, and is quite general for ordinary 
finite-size systems with such irrelevant interactions. This 
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property will be essentially changed by the introduction 
of non-hermitian terms, as is discussed in the next sec- 
tion. 

Before concluding this section, we wish to give an in- 
terpretation of the above spectral flow in terms of the 
Bethe-ansatz language At $ = the N rapidities 

for the ground state are distributed symmetrically on the 
real axis. As $ approaches 2ir, the last rapidity Am runs 
to infinity, while the others shift a little bit to the posi- 
tive direction, and are arranged in a configuration of the 
lowest energy state for the S z — 1/2 sector. This can be 



seen clearly if we divide the Bethe ansatz equations (2.6) 
at $ = 2ir into two classes, 



M-l 

(2.10) 

,M-1 and 

M-l 



for j = !,'■- 

Np(X M ) 



Y,6(\ M -\k) = 2ir{I M + l), (2.11) 



fc=i 



The last equation is satisfied by Am 
former equations reduce to 



+oo, and the 



M-l 

E 



Np(\j)+ V 6{Xj - A fc ) = 2*1$, 



where 



1 ' " 4 ' " ' ' 4 ' 



(2.12) 



(2.13) 



Beyond <& = 27r, the rapidity Am jumps from infinity 
to minus infinity and at last all rapidities return to the 
original configuration at $ = 4-7T. In this process, the 
rapidity Am crosses the branch cut of the logarithmic 
function in eq.(2.7), so that 



p(X^)=p(X^) + 2n. 



1. 



,M- 1, 



(2.14) 



This gives an interpretation of the period 47r, which is due 
to irrelevant interactions, in terms of the motion of the 
rapidities. In the next section, by extending the above 
idea, we discuss how the spectral flow changes its char- 
acter if we switch on the imaginary twist 



III. SPECTRAL FLOW FOR THE 
NON-HERMITIAN MODEL 

A. Basic properties of the spectral flow 



on the complex plane. In order to see clearly how the 
spectral flow occurs, let us regard the flow as a function 
of 4> with a fixed For small 'J, the flow is essentially 
the same as the $ = case, as can be seen in Fig. ||, 
where only the real part of E's is plotted for simplic- 
ity. Namely, except for the spectrum being complex, the 
qualitative nature of the whole spectral flow is not mod- 
ified by a finite but sufficiently small This implies 
that such small \& cannot collapse the l/N gaps due to 
irrelevant backward scatterings. It should be noted that 
at <f> = 2ir (mod 2ir) all eigenvalues are real. If \P is fur- 
ther increased, the gap between the first and the second 
excited state becomes small and eventually disappears at 
a certain \& = ^ CI , above which these two states become 
degenerate in the real part of E (They are complex con- 
jugate each other). Then the period of the ground state 
changes, as seen in Fig. |. In this case, the period of 
the ground state is 87r. If we further increase the pe- 
riod becomes longer. For example, in Fig. |[ it becomes 
macroscopic value 2ttN {= 12ir). Physically, the change 
of the period means that the small 1/iV gaps due to ir- 
relevant interactions are successively closed by the non- 
hermitian effects, as mentioned above. This phenomenon 
has the essentially same origin as for the metal-insulator 
transition triggered by non-hermitian hoppings. 




FIG. 1. Spectral flow of the S z = sector of the N = 6 with 
9 = system obtained by the exact numerical diagonalization 



1*3. 
OS 




FIG. 2. The same spectral flow but with * = 1.1 



We now extend the twist angle to a complex one. We 
should hence consider the spectral flow on various paths 



Let us now consider the spectral flow defined on the 
whole complex plane, which is sketched qualitatively 
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0.5 1 1.5 2 

FIG. 3. The same spectral flow but with * = 1.2 




0.5 1 1.5 2 

<Z>/2n 

FIG. 4. The same spectral flow but with * = 3 



in Fig. ![ If we start with the ground state at 6 = 0, 
and trace a path in plane, we may have such spectral 
flows as depicted by several curves in this figure. For 
small $ the path can be deformed to the real axis with- 
out changing the period. There exists, however, certain 
critical \l/ cr beyond which the gap between the first and 
the second excited states disappear. We denote it as a 
straight line in the figure. If we follow the path across 
this line, the period changes. It should be noted that the 
structure of Fig. || is periodic in $ = 2ir, though we have 
drawn only the lines relevant to the flow of the ground 
state starting from 9 = 0. Therefore, paths on the © 
plane are divided into several classes classified in terms 
of the periods of the ground state. Although the spectral 
flow on the complex plane has quite rich structure, their 
essential property is understood by studying the typical 
case where the period of the ground state changes from 
4tt to 8ir. In the following section, we will estimate the 
corresponding critical twist 'J'cr based on the finite-size 
corrections. 

Let us now check explicitly how the corresponding ra- 
pidities flow. The first example is Fig. [| which cor- 
responds to Fig. 0. Since we are concerned with 
the N = 6 site system, the ground state is described by 
three rapidities sitting near the real axis at $ = 0, indi- 
cated by dots. They continuously move on the complex 
plane, as <!> increases, towards the direction indicated by 
arrows. Finally at <f> — 2tt, rapidities reach the points 



2k 4n 6n 8% O 

FIG. 5. Schematic illustration of the spectral flow of the 
ground state starting at = 0. The ground state follows the 
curves towards the direction denoted by arrow and returns 
to the same ground state for the first time at the end of the 
curves. 
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FIG. 6. Rapidity flow of the ground state of the S z = 
sector with = 1.1, corresponding to Fig. |^. Initial ($ = 0) 
and final ($ = 2n) points are denoted by dots, flowing towards 
the direction indicated by arrows. 

indicated by dots. What should be noted here is the mo- 
tion of the last rapidity: After a long trip, it eventually 
reaches the imaginary axis. As a result, the whole dis- 
tribution of the rapidities is changed considerably, but 
it is still symmetric with respect to the imaginary axis, 
which ensures that the total energy is real. Contrary to 
this, for \P = 1.2 > ^ cr in Fig. |t] the final distribution is 
not symmetric, which suggests that there exists another 
state with the complex-conjugate energy corresponding 
to the configuration reflected with respect to the imagi- 
nary axis. We have checked that this is indeed true even 
for a larger system in Figs. ^ and 9. Therefore, we can 
see that whether the symmetric solution exists at $ = 2tt 
provides a criterion whether the gap exits. We have cal- 
culated for several systems, summarized in Table |. 

By using the criterion mentioned above, we wish to 
propose a practical method to estimate the critical value 
^ cr for a large but finite-size system. To this end let 
us assume that Am sits on the imaginary axis. Then if 
^ < ^cr, the Bethc ansatz equations on this assumption 
have a solution, while for '5 > fy CI , they do not have a 
solution. Therefore, the critical value of \& beyond which 
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FIG. 7. The same rapidity flow but with = 1.2, corre- 
sponding to Fig. |3| 



2.5, 

2 

1.5- 

1 

0.5- 





-0.5 0.5 1 1.5 
Re X 

FIG. 8. The same rapidity flow but of N = 20 system with 
* = 0.7 



there exist such solutions can be identified with ^ CI . The 
problem becomes simpler if we put it in a different way: 
We have so far solved the Bethe ansatz equations for 
a given ^. We now set the problem for seeking for a 
possible ^ provided that the value of imaginary Am is 
initially given. To solve this new problem, let us rewrite 
the Bethe ansatz equations at $ = 2ir, 



= 2-k% + l)+i*(A), 



(3.1) 
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FIG. 9. The same rapidity flow but of N = 20 system with 
* = 0.8 
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1.622 


1.759 
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^max 


1.0480 


0.9532 


0.8898 


0.7363 


0.6685 


0.5991 


A 


1.372 


1.446 


1.507 


1.712 


1.837 


1.999 


'I'jV.max 


1.1806 


1.0536 


0.9726 


0.7861 


0.7074 


0.6287 


A 


1.323 


1.398 


1.461 


1.673 


1.802 


1.967 



TABLE I. The values $ cr calculated by several ways. 
"V&cr" is the exact value calculated by the Bethe ansatz equa- 
tions. "$ mal " is the maximum v alue o f the function \I/(A) in 
(EO) by using the solution of Eq.fl2.12j). See III. B. "*jv, ma x" 



is the maximum value of the function ^jv(A) in (4.23) 



A/-1 

i#(A) = Np(iX) + e ( iX - x k) - 2tt(Im + 1) (3.2) 
fc=i 

In the above equations we should determine Xj (j = 
1, • • • , M — 1) and ^ for a given real A. The equations do 
not have a solution for $ > ^ cr so that for any given A, 
^(A) should not exceed ^cr. Namely, the function *(A) 
should have a maximum value, which in turn gives ^cr- 
Therefore, the problem is now to determine the function 
"J (A), and find its maximum, which indeed allows us to 
obtain the critical twist ^f r r- 



B. Bethe ansatz equations in a decoupling 
approximation 

Even if we employ the above trick, it is not so easy to 
determine ^ C r- We here propose an approximate treat- 
ment to the problem and confirm that it works fairly 
well for N = 20 system. Su ppose that N and A are 
large. Then the solutions of eq.( 2.12 ) can be regarded as 
approximate solutions. Therefore, if we substitute them 
into eq.(3.2) and plot the result as a function of A, extrap- 
olating to the small A region, we have, for example, Fig. 
|To| . This function indeed has a maximum value, which 
explains the existence of ^ CI . It is remarkable that the 
maximum value thus obtained is quite close to the ex- 
act one in Table |[ We list the maximum values of \17 (A) 
and corresponding A in the same Table. They deviate 
from the exact value by only a few percent in general. 
Therefore, we can say that the value of is essentially 
determined by the rapidity Aa/ on the imaginary axis, 
and our approximate treatment works quite well. 

In order to calculate ^> CT self-consistently, we have to 



calculate the finite-size corrections based on eqs.(3.1) and 
fl3.2p . There are various 0(l/N) terms, all of which 
should be taken into account. As we have already checked 
for the N = 20 system, however, the contributions from 
the third term of the left-hand side and the second term 
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A. Bulk limit 



FIG. 10. Function ff (A) in eq.(|3.2|) for N = 20 system. We 
use the solutions of eq.(2.12) as the nine rapidities A^ in this 
equation. 



i.e., 9(iX — Afc) and 
Therefore, we set 



of the right hand side in eq.(3.1 
ity(X), respectively, are quite small 
6{iX — Afc) — > 7r and i^(X) — » in this equation, and 
confirm the validity of this treatment by comparing the 
obtained results on this assumption with the exact ty cr 
for a small system. Thus, the Bethe equations for Xj 
(j = 1, • • • , M — 1) are decoupled from that for iX = Am, 
so that we use the notations fij for these Aj's. In order 
to proceed the calc ulatio n, let us introduce the following 
functions from eqs. (§T|) and @ 



Z N {n) 



pip) 



2ir 2nN 



N 



iV N (X) = Np(iX) + S{iX - /i fe ) - 7T + I J , (3.3) 

where Zjv(/Xfc) = I'k/N. All /i's are on the real axis. The 
distribution function is given by 



a Nil*) 



By definition, 



pM 

2tt 



2ttA 



d/j.cTN(p) = 



1 

AT 



(3.4) 



(3.5) 



The above equations enable us to further proceed the 
analytic calculation, which will be performed separately 
in the next section. 



IV. FINITE-SIZE CORRECTIONS 

We are now ready to evaluate finite-size corrections 
based on the equations in the preceding section. In par- 
ticular, we derive the ^ CT as a maximum value of the 
function W(A) in the bulk limit with leading finite-size 
corrections. So far, practical methods to compute the 
finite-size corrections for non-hermitian models [ p8|j30| 
have not been well explored. So, we hope that the follow- 
ing treatment may also shed some light on such problems. 



We wish to first investigate Eqs. (3. 4) and (3.3) in the 
bulk limit N — > oo to confirm that our approach is indeed 
consistent. Since the effects of the twist angle should be 
0(1/N), we expect that the function ^(A) in the bulk 
limit is identically zero for any A, which is shown to 
be in dee d satis fied in our treatment. In the bulk limit, 
Eqs.(O) and Q 



Coo(A0 = 
^oo(A) 

N 



P'jp) 
= p(iX) 



are 
1 

27 



dfj,'9'(fx - //)o-oo(m' 

oo 

d/i9(iX — jj)a 



(4.1) 
(4.2) 



where the left-hand side of the second equation means 



limjv- 



N 



Since /i's are distributed on the real 



axis, their distribution can be easily determined by the 
Fourier transformation. The result is 

Coo(m) = ^scch7r^. (4.3) 

By substituting this solution into eq.( [4.2| ), we have 

i^oo(X) [°° , sech7ru ( li — i(X— 1 
' ilji In ' 



N 



-i In 



A 



2 
1/2 



fi-i(X + l) 



A + 1/2 



(4.4) 



where we have assumed A > 1. The first term of the right- 
hand side is evaluated by extending the contour integral 
in the lower half plane and summing up the infinite scries 
of residues on the imaginary axis. We find that it cancels 
the second term, and leads to 



i^oo(A) 
N 



= 0. 



(4.5) 



This proves that our treatment indeed produces the cor- 
rect result in the bulk limit because effects of the twist 
should disappear in this limit. 



B. 1/N corrections 

So far we have checked that ^oo(A) = for any A > 1, 
which means that the gap in the spectrum always disap- 
pears in the bulk limit. It is quite natural because the 
gap is caused by marginally irrelevant interaction. Now 
we proceed to finite-size systems, for which we expect a 
finite ^cr. 

Subtracting eq. 



Jj) from eq.(3.4) we have the identity 
valid for a system with arbitrary N, 

/>OG 

crjv(iu) - <Too(/i) = - / d/j,'q(n - n') 



, (4.6) 
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where the kernel is 
20) 



due- luJfl - 



2tt ./_ 0O eM + 1 



(4.7) 



Denote the largest rapidity on the real axis as A. Then 
Ztv(A) = i — jf, or equivalently, 



d^tJAr(/i) 

Similarly, \P(A) is rewritten as 



1 2 

2 ~ TV' 



(4.8) 



AT 



dfiQ(\,(j,) 



(4 ' J) 



where 



Q(A, ^) = 0{i\ -n)- / d/*'0(*A - f/)q(p' - »), 



(4.10) 

The kernel <3(A, /1) can be evaluated by counter integral 
in the lower half plain, yielding the result 



Q(A,/i) = -- 
Q(X,fi) = iln 



-Q(A,/i), 

fi-i(\- 1) 
/i — zA 



(4.11) 



Substituting this equation into Eq.(4.9), we find that the 
first constant term in Eq.(4.11) cancels the last constant 
term in Eq.(4.9). Therefore, 



^jy(A) 
N 



dfj,Q(X,fx) 

o 



(4.12) 



To evaluate the leading finite-size effects, let us apply 
the Euler-Maclaurin formula to Eqs.(4.6) and (4.9). Up 
to 0(N~ 2 ), we have 



d^'q{n - n')a N {n') 



q{^-A) q'ifx-A) 



2N 



12N 2 a N (A) 



(-A terms), (4.13) 



N ja 

, Q(A,A) , Q'(A,A) 



N 



6N 2 a N (A)' 



(4.14) 



where Q'(X, A) = dQ(X,A)/dA. The cut-off A is deter- 
mined by the condition 



/•OO g 



(4.15) 



_Z. Finite size corrections to aff(n) 

As we have exploited the decoupling approximation, 
we can evaluate the finite-size corrections to <7jv(A) fol- 
lowing standard techniques |3l|-|34[| . Introducing 



o-jv(m) = 



for fj, < 0, 

for /u > 

cjvCj" + A) for /u < 0, 



we can write down eq.( 4.13| ) in the region fi 
neglecting —A terms, 



df/q(jjt- fj,')a%(jjt') 



(4.16) 
A, by 



a+(fi)+a N {fi) - 



2N 12Af 2 CTAr(A) ' 



(4.17) 



Solving the equation in a standard way summarized 
briefly in the appendix, we end up with the final results 

a+(co) = C(uj) + G+{u/2k) [P(uj) + Q+(w)] , (4.18) 

where a^(u>) is the Fourier-transformation of cr^(A) and 
where 



C[w) 



1 



1 



2n\2n 12AT 2 aiv(A)A 



G + (cj/27t) = V2W 

' 2ne 



/r ( i — — 

7 1 2 2tt 



P( W ) = 
with 



1 



1 



2nig 



/2^\2N 12N 2 a N {A) 12N 2 o N {A) J 

-ttA 



1 e 



'2n w + 7rz 



*5 = 1/24, 



7 G+ (i/2) . 



(4.19) 



ojv(A) = —a, 



A/ 



A = - In 

7T 



7re/3 



(4.20) 



Here several constants are given by 



1 



a = - ( a + \/al - 4ai 



= 



a = 



«i 



1 1 



2 V2 6a 

1 1 

2 + 7f 4 ' 

^ + (ig) 2 



ig, 



(4.21) 
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2. Finite-size corrections to ^jv(A) 

By using the preceding results for tJAr(A), let us now 
calculate ^/at(A) with finite-size effects. To begin with, 
we conv ert the first term in the right-hand side of 
eq.( 4.14 ) into the integral in Fourier space, 



-2 / dnQ{\,/j,)a N (ii) 
I A 



'2tt 



due 



-AM 



-(A-l)M 



a?+(w). (4.22) 



It is easily shown that C(lo) in a^(u>) cancels o ut th e 
second and third terms in the right-hand side of eq.( 4.14 ). 
Therefore, we eventually arrive at 



*iv(A) 



-(A-l)M 



due 



xG+{uj/2it)(P(lu) + Q+(lu)). 



(4.23) 



We have evaluated the integral numerically for the 
N = 20 system in Fig. O, corresponding to Fig. [!(]. 
As we expect, the leading finite-size calculation yields 
the desirable results. In fact, the maximum value of this 



0. 7 



0.6 



0.5 



).4 




FIG. 11. *jv(A) for N = 20 system 

function reproduces the critical \P expected for small sys- 
tems fairly well, which are listed in Table Q. 




200 



FIG. 12. The data of in Table || and their least-square 
fits. See the text. 

It seems not easy to precisely determine the large- 
N behavior of the function ^E'(A). However, observ- 
ing the data for several small N systems of the exact 



Bethe ansatz calculation and also for larger N systems of 
Eq.(4.23), we can predict 1/ In N behavior in the asymp- 
totic regime. To confirm this explicitly, we plot in Fig. |l2| 
the least-square fit using (In TV) - ™ with n = 1,2,3 with 
the data points. The solid- line is for the exact Bethe 
ansatz results, while the dotted- and dashed-lines are for 
the decoupled Bethe ansatz results of finite-size systems 



in section III. B. and for Eq.(4.23), respectively. As ex- 
pected, the larger the system is, the smaller the difference 
between the dotted- and dashed-lines becomes. Further- 
more, these \l/ cr in a decoupling approximation are very 
close to the exact \l/ cr . Note that \& cr is not a gap itself 
between the first and second excited states. The gap of 
0(1/ N) exist there, and we need ^ cr of 0(1/ In N) to 
close this gap. 



V. SUMMARY AND DISCUSSIONS 

Stimulated by the recent advance in the theory 
for localization-delocalization transitions in ID non- 
hermitian systems, we have studied how the non- 
hermitian property affects the integrable Heisenberg spin 
model. In particular, we have focused the particular at- 
tention on the spectral flow by imposing twisted bound- 
ary conditions with complex twist. By solving the Bethe 
ansatz equations numerically for finite systems, we have 
shown that a global structure of the spectral flow com- 
pletely changes its character as the imaginary twist in- 
creases. Namely, the period of the ground state increases 
from 47T to 8ir at a certain imaginary twist angle. If 
the imaginary twist is further increased, the period suc- 
cessively jumps to a higher value at the corresponding 
imaginary twist angle. The origin of this phenomenon 
has been clarified: the small 1/N energy gaps for a fi- 
nite system caused by irrelevant interactions are closed 
by the non-hermitian hopping of XY term. To determine 
the critical imaginary twist for the transition, we have 
also proposed a practical method to evaluate the finite- 
size corrections, and discussed the behavior in the large 
N regime. Although we have neglected some terms which 
may give 0(1/N) contributions in this treatment, the es- 
timated value of 'J'cr is quite close to the exact one. Since 
practical ways to explicitly evaluate the finite-size correc- 
tions for non-hermitian models have not been developed 
well, the present treatment may provide a guideline to 
study such problems in more detail. 

In this paper we have treated $ and \& on an equal foot- 
ing. Before concluding the paper, we briefly make a com- 
ment on the physical origin of these twists. As mentioned 
in the introduction, the real twist 4> can be naturally re- 
garded as an Aharonov-Bohm flux. As for imaginary \E', 
its origin may be ascribed to randomness in the system. 
This is indeed seen by observing the Hamiltonian, 
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H 



E 



i c+ c- 



with the boundary condition 



S 



j+N 



o QZ nz 

J ji '-'j+N — D j ■ 



(5.1) 



(5.2) 



where 9j are not necessarily uniform. The Bethe ansatz 
equation is the same as Eq.([2.6|), but with 



B = <P + iJ2^i- 



(5.3) 



the mesoscopic gaps, whereas that for the XXZ case is 
by the collapse of the ordinary gap. We can now see that 
there are deep relationship among the transition found 
in the spectral flow of the present model, the "insulator- 
metal" transition for the XXZ model, and the Hatano- 
Nelson model, although at first sight our results may not 
have any essential connections to latter phenomena; they 
may indeed share some essential features, e.g., the transi- 
tion is driven by imaginary vector potentials and is char- 
acterized by complex eigenvalues at criticality. It may 
be interesting to study wider class of models with imagi- 
nary vector potentials to understand a novel class of the 
metal-insulator transition in a unified framework. 



Therefore, it is seen that our imaginary twist can be re- 
garded as the consequence of non-hermitian random hop- 
pings. By taking into account the ensemble average over 
random variables , we may have a chance to discuss the 
effects of random non-hermitian hoppings, which should 
be treated in the future work. 



The Hamiltonian (5.1) is also useful to see how the 
present analysis is related to the metal-insulator tran- 
sition in the asymmetric models, and also to the der- 
ealization transition in the Hatano-Nelson model. To 
this end, let us set = *f?/N, where vff is given by 



(2.3). By using the hard-core boson representation, the 
XY term is converted into the hopping term of bosons, 



#hop = £,(e-*/ w 6W + e */"&t 



where bj is a 



3 y 3 

boson operator. Namely, the present model is equivalent 
to the interacting bosonic lattice model with an imagi- 
nary vector potential. Note that the kinetic part is es- 
sentially the same as that for the Hatano-Nelson model. 
For the latter model, Hatano and Nelson demonstrated 
that the "insulating" state caused by randomness can be 
delocalized and become metallic at a certain critical '3/ 
when the strength of the imaginary vector potential is 
increased. This implies that the transition is caused by 
the competition between the imaginary vector potential 
and the random potential. We can also observe a simi- 
lar insulator-metal transition caused by imaginary vector 
potential for the periodic model such as the asymmetric 
XXZ model with Ising-type anisotropy , for which the 
massive (insulating) state is formed by the many-body 
interaction, instead of randomness. In contrast to the 



above examples, for the isotropic Heisenberg model ( |5.l| ) 
we have dealt with, the system is massless and "metallic" 
in the thermodynamic limit even for $ = 0. However, 
we should recall that there still exist small (mesoscopic) 
gaps formed by irrelevant many-body interactions for a 
finite-size system, which control the period of the spec- 
tral flow for the ground state. We have shown explicitly 
via the present study that the competition between the 
imaginary vector potential and the many-body interac- 
tion indeed modifies the period, which implies that the 
mesoscopic gaps are closed by the imaginary vector po- 
tential: Our transition is characterized by the collapse of 
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APPENDIX A: WIENER-HOPF EQUATION 

Since there exist extensive studies of finite-size calcula- 
tions we summarize it briefly for reference sake. 

Introducing the Fourier transformation with the nor- 
malization 



4(A) 



1 



duje- lXuJ a±(uj), (Al) 



Eq.(4.17) is given by the Fourier component 

I — 

2 C(w) (A2) 



NK " J ' e-H + l 2V2n 

The kernel is now factorized as 

e - H + l = G + (w/27r)G_(w/27r), (A3) 

where G±(z) is holomorphic in the upper-(lower) half 
plain, respectively. The asymptotic form is 



z 2z A 

Then Eq. ([A2]) can be written as 



(A4) 



-Q+H 



G+(lo /2tt) 
= -G-{u/2n){a N {iu) + C{w)) + Q_(w) 
= P(u). (A5) 
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where Q±(u>) is holomorphic in the upper-(lower) half 
plain, satisfying 



G-{uj/2-k)- 



v sech£ 



2V2vr 



= Q+H + Q_H, (A6) 



where we take into account only the contribution from 
the pole nearest to the real axis in the text. From the 
asymptotic form Eq.(A4), P(ui) is determined, which is 
summarized in the text. 

The parameters A and c(A) are determined by the 
conditions Eq.(4.15) and 



<t n {A) = 



or more precisely, 
G+(z/2) e ^ A : 



dwcr^(w), 



7r n{2iTig) 



(A7) 



2 s/2) N l2N 2 a N {A)' 
{2mg) (2irig) 2 



2N 2AN 2 a N (A)' 
The solutions are summarized in the text. 
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